Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 13 December 2012 (MN MbX style file v2.2) 



Measurements of H{z) and Da{z) from the Two-Dimensional 
Two-Point Correlation Function of Sloan Digital Sky Survey 
Luminous Red Galaxies 

Chia-Hsun Chuang* and Yun Wang 

, Homer L. Dodge Department of Physics & Astronomy, Univ. of Oklahoma, 440 W Brooks St., Norman, OK 73019, U.S.A. 



13 December 2012 



ABSTRACT 

We present a method for measuring the Hubble parameter, H{z), and angular diameter 
distance, Da{z), from the two-dimensional two-point correlation function, and validate it us- 
ing LasDamas mock galaxy catalogs. Applying our method to the sample of luminous red 
galaxies (LRGs) from the Sloan Digital Sky Survey (SDSS) Data Release 7 (DR7), we mea- 
sure H{z = 0.35) = ff(0.35) = 82.11^9 kms-iMpc"\ Da{z = 0.35) = L>a(0.35) = 
lOAStfs Mpc without assuming a dark energy model or a flat Universe. We find that the de- 
rived measurements of H{0.35) rs{zd)/c and DA{Q-i5)/rs{zd) (where rs{zd) is the sound 
horizon at the drag epoch) are nearly uncorrected, have tighter constraints and are more 
robust with respect to possible systematic effects. Our galaxy clustering measurements of 
{H{Q.35) rs{zd)/c, DA{0.35)/rs{zd)} = {0.0434 ± 0.0018, 6.60 ± 0.26} (with the corre- 
lation coefficient r = 0.0604) can be used to combine with cosmic microwave background 
and any other cosmological data sets to constrain dark energy. Our results represent the first 
measurements of H{z) and Da{z) (or H{z)rs{zd)/c and I?^(0.35)/rs(2:d)) from galaxy 
clustering data. Our work has significant implications for future surveys in establishing the 
feasibility of measuring both H{z) and Da{z) from galaxy clustering data. 
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1 INTRODUCTION 

The cosmic large-scale structure from galaxy redshift surveys pro- 
vides a powerful probe of dark energy and the cosmological model 
that is h ighly complementary to the cosmic micr owave background 
(CMB) /Benn ett et ai]|2003h . supemovae (SNe) ("Riess et al 199^; 



Perlmutter et alT 199S ), and weak lensing (Wittman et al. 200| 



Bacon, Refregier, & Ellis' I2000I; iKaiser. Wilson, & Luppinol 12001 



IOCS; 



Ivan Waerbeke et al.. .2000i) . The scope of galaxy redshift surveys 
has dramatically increased in the last decade. The PSCz surveyed ~ 
15, O OP galaxies using the Infrared Astronomical Satellite (IRAS) 
([senders et al. 2000), the 2dF Galax y Redshift Survey (2dF GRS) 
obtained 221,414 galaxy redshifts ( Colless et al1l200Ul2003h . and 
the Sloan Digital Sky Survey (SDSS) has c ollected 930,000 galax y 
spectra in the Seventh Data Release (DRV) jAbazaiian et alj2009h . 
The ongoing galaxy surveys will probe the Universe at higher 
redshifts; WiggleZ is surveying 240,000 e mission-line galax ies at 
0.5 < z < 1 over 1000 square degrees jBlake et al.ll2009l) , and 
BOSS is surveying 1.5 million luminous red galaxies (LRGs) at 
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0.1 < 2 < 0.7 over 10,000 square degrees dEisenstein et al.l201 ih . 
The planned space mission Euclid will survey over 60 million 
emission-line galax i es at 0.5 < z <2 over 20,000 square degrees 
dCimatti et al.ll2009l ; IWang et al.ll2010h . 

Large-scale structure data from galaxy surveys can be an- 
alyzed using either the power spectrum or the correlation func- 
tion. Although these two methods are simple Fourier trans- 
forms of one another, the analysis processes are quite differ- 
ent and the results cannot be converted using Fourier trans- 
form directly because of the finite size of the survey vol- 
ume. The SDSS data have been analyzed using bot h the 
powe r spectrum method (see, e.g., Tegmark et al. 2004^ 



Hutsi 



2005; Padmanabhan et al. 2007; Blake et ali 120071 : iPercival et al 
2007, 2010; Reidetal. 2009), and the correlation function 



method (see, e.g., feisenstein et al.l l2005l: lOkumura et al.l 20081: 
Cabre & Gaztanaga' '2008'; 'Martinez et al' '2009"; 'Sanchez et alj 
2009; Kazinetal. 2010a; Chuang, Wang, & Hemantha 201$^ 
While previous work has focused on the spherically averaged two- 
point correlation function (2PCF), or the radial projection of the 
two-dimensional two point correlation function (2D 2PCF), we 
measure and analyze the full 2D 2PCF of SDSS LRGs in this study. 
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The power of galaxy clustering as a dark energy probe lies 
in the fact that the Hubble parameter, H{z), and the angular di- 
ameter distance, Da{z), can in principle be extracted simultane- 
ously from data through the measurement of the baryon acous- 
tic oscillat ion (BAO) scale in the radial and transverse direction s 
( iBlake & Glazebrook 20011 : ISeo & Eisensteinll2003l : IWang|l2006h . 
This has not been achieved in th e previous work in the analysis of 
real data. lOkumura etlTI (|2008[) concluded that SDSS DR3 LRG 
data were not sufficient for measuring H{z) and Da{z); they de- 
rived constraints on cosmologica l parameters assuming tha t dark 
energy is a cosmological constant. ICabre & Gaztanagal (l2008h mea- 
sured the linear redshift space di stortion parameter /3, galaxy bias , 
and CT8 from SDSS DR6 LRGs. lOaztanaga. Cabre. & Huil j2009h 
obtained a measurement of H{z) by measuring the pea k of the 
2PCF along the line of sight. However. lKazin et aL H2010hl) showed 
that the amplitude of the line-of-sight peak is consistent with sam- 
ple variance. 

In our previous paper, *Chuang. W ang. & Hemanthal ( l2012h . 
we presented the method to obtain dark energy and cosmologi- 
cal model constraints from the spherically-averaged 2PCF, with- 
out assuming a dark energy model or a flat Universe. We demon- 
strated the feasibility of extracting H{z) and Da{z) by scaling 
the spherically-averaged 2PCF (which leads to highly correlated 
measurements). In this paper, we obtain robust measurements of 
H{z) and Da{z) through scaling, using the 2D correlation func- 
tion m easured from a sample of SDSS DR7 LRGs jEisenstein et al.l 
l200lh . This sample is homogeneous and has the largest effec- 
tive survey vol ume to date for studying the quasi-linear regime 
jEisenstein et aL 2005,) . In Section|2l we introduce the galaxy sam- 
ple used in our study. In Section [3] we describe the details of our 
method. In Section|4l we present our results. In Sectioin|5] we ap- 
ply some systematic tests to our measurements. We summarize and 
conclude in Sec. [6] 



2 DATA 

The SDSS has observed one-quarter of the entire sky and per- 
formed a redshift survey of galaxies, quasars an d stars in five 
passbands w, g, r, i, and z w ith a 2.5m telescope jpukugita et al.l 
ll996l : lGunn et al.ll998ll20o3) . We use the public catalog, the NYU 
Value-Added Galaxy Catalog (VAGC) (Blanton et al. 2005), de- 
rived fro m the SDSS II final public data release. Data Release 7 
(DR7) tAbazaiian et all 12009 ). We select our LRG sample from 
the NYU VAGC with the flag primTarget bit mask set to 32. 
K-corrections have been applied to the galaxies with a fiducial 
model (ACDM with flm = 0.3 and h = 1), and the selected 
galaxies are required to have rest-frame g-band absol ute magni- 
tudes -23.2 < Mg < -21.2 telanton & Roweij2007h. The same 
selection criteria were used in previous paper s (IZehavi et al ]|2005l : 
lEisenstein et"ai] |2005| ; [Okumura et al] l2008l : |Kazin et al." 2010a'). 
The sa mple we use is referred to as "DR7fuir' in Kazin et al. 
( l2010ij) . Our sample includes 87000 LRGs in the redshift range 
0.16-0.44.0 

Spectra cannot be obtained for objects closer than 55 arcsec 
within a single spectroscopic tile due to the finite size of the fibers. 



To coiTect for these "collisions", the redshift of an object that failed 
to be measured would be assigned to be the same as the nearest 
successfully observed one. Both fiber collisi on corrections and K - 
corrections have been made in NYU- VAGC telanton et alll2005h . 
The collision corr ections applie d here are different from what has 
been suggested in lZehavi et al.l 12005) . However, the effect should 
be small since we are using relatively large scale which are less 
affected by the collision corrections. 

We construct the radial selection function as a cubic spline fit 
to the observed number density histogram with the width Az = 
0.01. The NYU- VAGC provides the description of the geometry 
and completeness of the survey in terms of spherical polygons. We 
adopt it as the angular selection function of our sample. We drop 
the regi ons with co mpleteness below 60% to avoid unobserved 
plates teehavi et al.ll2005h . The Southern Galactic Cap region is 
also dropped. 



3 METHODOLOGY 

In this section, we describe the measurement of the correlation 
function from the observational data, construction of the theoret- 
ical prediction, and the likelihood analysis that leads to constraints 
on dark energy and cosmological parameters. 



3.1 Measuring the Two-Dimensional Two-Point Correlation 
Function 

We convert the measured redshifts of galaxies to comoving dis- 
tances by assuming a fiducial model, ACDM with fl™ = 0.25. 
We use the two-point correlation function estimator given by 
lLandv&Szaiavl ( ll993h : 



DD{a, n) - 2DR{g, tt) + RR{g, tt) 



(1) 



where tt is the separation along the light of sight (LOS), a is the 
separation in the plane of the sky, DD, DR, and RR represent 
the normalized data-data, data-random, and random-random pair 
counts respectively in a distance range. The LOS is defined as the 
direction from the observer to the center of a pair. The bin size 
we use in this study is 10 ^~^Mpcx 10 /i~^Mpc. The Landy and 
Szalay estimator has minimal variance for a Poisson process. Ran- 
dom data are generated with the same radial and angular selection 
functions as the real data. One can reduce the shot noise due to ran- 
dom data by increasing the number of random data. The number of 
random data we use is 10 times that of the real data. While calculat- 
ing the pair counts, we assign to each data point a radial weight of 
1/[1 + n{z) ■ Ptu], wher e n(z) is the radial sele ction function and 
P„ = 4 • 10* /i"^Mpc^ teisenstein e7Zll2005h . We use the same 
as Eisenstein et al. (2005) in which they used the sample of the 
SDSS DR3. Although the data release versions are different, the 
properties of the galaxy sample should be basically the same. We 
find that the error bars estimated from LasDamas mock catalogs 
could be improved by 10% while using the weighting (compared 
to the error bars obtained without using the weighting). We expect 
that the results should not be sensitive to the used. 



^ We have identified a bug while computing the weighting of each galaxy 
in the first draft of this paper. The bug was that we computed the weights of 
random data with the number density of random data instead of observed 
data. Tills would introduce a bias when the number density is not homoge- 
neous. 



3.2 Theoretical Two-Dimensional Two-Point Correlation 
Function 

We co mpute the linear power spectra at z — 0.35 by using 
CAMB jLewis. Challinor. & Lasenbvll2OO0l) . To include the effect 
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of non-linear structure formation on the BAOs, we first calculate 
the dewiggled power spectrum 



Pdw{k) = P!i„(fc)exp ( )+Pn«,(fe) 



exp 



,(2) 



where Pun{k) is the linear matter power spectrum, P„,„ (fc) is the 
no-w iggle or pure CDM pow er spectrum calculated using Eq.(29) 
from lEisenstein & Hul l ll998l) . and fc* is marginalized ovei0 with a 

flat prior over the range of 0.09 to 0. 13 /iMpc~^^ 

We then use the software package halofit jSmith et"al]|2003l) 
to compute the non-linear matter power spectrum: 



rhaiofitik) 



Pjialo f it ,nw (^) 



(fc) 

Pnl{k) = Pdw{k)rhalofit{k), 



(3) 
(4) 



where Phaiofit,nw{k) is the power spectrum obtained by applying 
halofit to the no-wiggle power spectrum, and Pni{k) is the non- 
linear power spectrum. We compute the theoretical real space two- 
point correlation function, ^(r), by Fourier transforming the non- 
linear power spectrum P^i (k) . 

In the linear regime (i.e., large scales) and adopting the small- 
angle approximation (which is valid on scales of interest), the 2D 
correlation function in the redshift space can be written as jKaised 
ll987l ; lHamiltoi]ll992l) 

r('T,^) =Co(s)Po(m)+6(s)P2(m)+C4(s)P4(^), (5) 

where s — \/a'^ + tv^, fi is the cosine of the angle between 
s = {a, tt) and the LOS, and Pi are Legendre polynomials. The 
multipoles of ^ are defined as 

Co(r) = 



35 



7 



where /3 is the redshift space distortion parameter and 



^(r')r"^dr' . 



(6) 
(7) 
(8) 

(9) 
(10) 



Next, we convolve the 2D correlation function with the distribution 
function of ra ndom pairwise velocities, f{v), to obtain the final 
model tt) l |Peebleslll980l) 



^(a,7r) = / ^ (T,7r - 



H{z)a{z] 



f(v)dv, 



(11) 



where the random motions are rep resented by an exponential form 
( iRatcliffe et al.lll998l : lLandvll2002h 



(Ti,\/2 



exp 



(12) 



^ Although kj, can be computed by renorma lization perturbation theory 
iCrocce & Scoccimarrol2006tlMatsubaral2007l) . doing so requires knowing 
the amplitude of the power spectrum, which is also marginalized over in 
this study. 



where a„ is the pairwise peculiar velocity dispersion. 

The parameter set we use to compute the theoretical cor- 
relation function is {H{z), Da{z), l3,Q,m,h? ,^bh^ ,ris,(Jv,ki,}, 
where Q,m and fi;, are the density fractions of matter and baryons. 
Us is the powerlaw index of the primordial matter power spec- 
trum, and h is the dimensionless Hubble constant {Ho = lOO/i 
km s~^Mpc~^). We set h — 0.7 while calculating the non-linear 
power spectra. On the scales we use for comparison with data, the 
theoretical correlation function only depends on cosmic curvature 
and dark energy through parameters H(z) and Da{z), assuming 
that dark energy perturbations are unimportant (valid in simplest 
dark energy models). Thus we are able to extract constraints from 
data that are independent of a dark energy model and cosmic cur- 
vature. 

Fig pTa)] shows the 2D 2PCF measured from SDSS LRGs 
compared with a theoretical model. The measured 2D 2PCF of the 
SDSS LRGs has been smoothed by a Gaussian filter with rms vari- 
ance of 2/i^^Mpc to illustrate the comparison of data with model 
in this figure, as the unsmoothed data are very noisy. Smoothing 
is not used in our likelihood analysis to avoid possibly introducing 
systematic biases. Fig. 1 1(b) [ shows the 2D 2PCF measured from a 
single LasDamas SDSS LRG mock catalog for comparison. The 
similarity between the data and the mock in the range of scales we 
used (indicated by the shaded disk) is apparent. 

Fig|2] shows the averaged 2D 2PCF measured from the Las- 
Damas mock catalogs compared with a theoretical model. The con- 
tour levels are apparent in the measured 2D 2PCF even though 
no smoothing is used; this is due to the reduction of shot noise 
achieved by averaging over 160 mock catalogs. Clearly, our 2D 
theoretical model provides a reasonable fit to data on intermediate 
(and quasi-linear) scales. The deviations on smaller scales may be 
due to the simplicity of the peculiar velocity model we have used. 
We do not use the smaller scales (s < 40 /i^^Mpc), where the scale 
dependence of redshift distortion and galaxy bias are not negligi- 
ble and cannot be accura t ely de termined at p resent. According t o 
Fig. 5 in'E isenstein et al.l ( l2005h and Fig. 4 in lBlake et all ( 1201 1[) . 
these effects are negligible at s > 40 fe~'^Mpc. On large scales, 
data become very noisy as sample variance dominates. For these 
reasons, we will only use the scale range of s = 40 — 120 /i^^Mpc 
in our analysis. We do not consider wide-angle effects, since 
they hav e been shown to be small on the length scales of inter- 
est here jSamushia et al. 1 1201 ih . ISam ushia et ah] ( 1201 ll) showed 
that the corrections (i.e. nonlinear effect and wide-angle effect) 
to the Kaiser formula are small comparing to the statistical er- 
rors on the measurement of the correlation function from SDSS 
DR7 LRG for the scale range interested (s=40-120 Mpc/h). In this 
study, we include the largest correction, dewiggling (nonlinear- 
BAO), and the nonlinear effects at small scales. Since including 
a larger range of scales gives more stringent constraints, our choice 
of s = 40 — 120 /i~^Mpc represents a conservative cut in data to 
reduce contamination by systematic uncertainties. 



3.3 Covariance Matrix 

We use the mock catalogs from the LasDamas simulation^ 
(McBride et al., in preparation) to estimate the covariance matrix 
of the observed coiTelation function. LasDamas provides mock cat- 
alogs matching SDSS main galaxy and LRG samples. We use the 
LRG mock catalogs from the LasDamas gamma release with the 
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(a) 2D 2PCF from SDSS LRGfuU sample 



(b) 2D 2PCF from single mock catalog 



Figure l. |(a)| The two-dimensional two-point correlation function (2D 2PCF) measured from SDSS DRV LRGs in a redshift range 0.16 < ^ < 0.44 (solid 
black contours), compared to a theoretical correlation function with parameters close to the best fit values in the likelihood analysis (dashed red contours). 
The theoretical model has H(z = 0.35) = SLSkms"! Mpc-\ Da{z = 0.35) = 1042 Mpc, (3 = 0.35, Qmh'^ = 0.117, ftth^ = 0.022, = 0.96, 
(jy = 300km s~^, and fe* = 0.11. |(b)| The 2D 2PCF measured from a single mock catalog, compared to a theoretical model with the input parameters of 
the LasDamas simulations and {/3, (t„, fc*} are set to {0.35, 300km s~^, 0.11/iMpc^^} (dashed red contours). In both figures, the shaded disk indicates the 
scale range considered (s = 40 — 120 ft,~^Mpc ) in this study. The thick dashed blue circle denotes the baryon acoustic oscillation scale. The observed 2D 
2PCF has been smoothed by a Gaussian filter with rms variance of 2/i~^Mpc for illustration in these figures only; smoothing is not used in our likelihood 
analysis. The contour levels are § = 0.5, 0.1, 0.025, 0.01, 0.005, 0. The § = contours are denoted with dotted lines for clarity. 
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Figure 2. The average two-dimensional two-point congelation function 
(2D 2PCF) measured from 160 LasDamas SDSS LRGfuU mock cata- 
logs (solid black contours), compared to a theoretical model with the in- 
put parameters of the LasDamas simulations and {/3, cr„, fcj,} are set to 
{0.35,300km s-\0.11/iMpc"i} (dashed red contours). The gray area 
is the scale range considered (s = 40 — 120 h~^Mpc ) in this study. 
The thick dashed blue circle denotes the baryon acoustic oscillation scale. 
The contour levels are apparent in the 2D 2PCF measured from mock 
catalogs, even though no smoothing is used. The contour levels are ^ = 
0.5, 0.1, 0.025, 0.01, 0.005, 0. The ^ = contours are denoted with dot- 
ted lines for clarity. 



same cuts as the SDSS LRG DR7full sample, -23.2 < A4g < 
—21.2 and 0.16 < z < 0.44. We have diluted the mock catalogs 
to match the radial selection function of the observed data by ran- 
domly selecting the mock galaxies according to the number density 



of the data sample. We calculate the 2D correlation functions of the 
mock catalogs and construct the covariance matrix as 

JV 

fc=i 

where is the number of the mock catalogs, is the mean of the 
m"' bin of the mock catalog correlation functions, and is the 
value in the m"* bin of the fc*'^ mock catalog correlation function. 
Note that the covariance matrix constructed from the LasDamas 
mock catalogs is noisy because only 160 mock catalogs are avail- 
able. Therefore, we smooth it before using it to calculate the likeli- 
hood. 

3.4 Likelihood 

The likelihood is taken to be proportional to exp(— x^/2) 
JPress et ^1 19921) . with given by 

where Ntins is the number of bins used, Sm = {<ym,Trm), 
(,th is the theoretical correlation function, and ^ots is the 
observed correlation function. Note that ^th(Si) depends on 
{H{z), Da{z),P, Q.mh'^, D,bh'^,ns,a^,k-t}. 

In principle, we should recalculate the observed correlation 
function while computing the for different models. However, 
since we don't consider the entire scale range of the correlation 
function (we only consider s = 40 — 120 h~^Mpc in this study), 
we might include or exclude different data pairs for different mod- 
els which would render values arbitrary. Therefore, instead of 
recalculating the observed correlation function, we apply the in- 
verse operation to the theoretical correlation function to move the 
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parameter dependence from the data to the model, thus preserving 
the number of galaxy pairs used in the likelihood analysis. 

Let us define T as the operator converting the measured cor- 
relation function from the fiducial model to another model, i.e., 



(15) 



where (^^^^ (s) is the observed correlation function assuming the 
fiducial model. This allows us to rewrite as 



E 



{r-M6/.(s.)]-d"fe)}, 



(16) 



where we have used Eqs.dlBt and J15t . 

To find the operator T, note that the fiducial model is only 
used in converting redshifts into distances for the galaxies in our 
data sample. In the analysis of galaxy clustering, we only need the 
separation of a galaxy pair, and not the absolute distances to the 
galaxies. For a thin redshift shell, the separation of a galaxy pair 
in the transverse direction is proportional to Da{z)A6 (A6 is the 
angle separation of the galaxy pair) and the separation along the 
direction of the line of sight is proportional to Az/H{z) (Az is the 
redshift difference between the galaxy pair). Thus, we can convert 
the separation of one pair of galaxies from the fid ucial model to an- 
other m odel by performing the scaling (see, e.g.. lSeo & EisensteinI 
(Eooi)) 



(a',7r') 



DAjz) IP 

D{:'{z)''' H{z) 



(17) 



Therefore, we can convert the measured 2D correlation function 
from some model to the fiducial model as follows: 



Da{z) H 
a 



fid 



Di^'iz) 



H{z 



(18) 



This mapping defines the operator T. 

We now apply the inverse operation to the theoretical correla- 
tion function: 



-i (c,th[(y,TV)) = t,th I — — cr, 



can be calculated by substituting eq. l |19l l into eq. l |16l >. 



(19) 



3.5 Markov Chain Monte-Carlo Likelihood Analysis 

We use CosmoMC in a Marko v Cha in Monte-Carlo 
likelihood analysis l lLewis & Bridi3 l2002h . The param- 
eter space that we explore spans the parameter set of 



{ii"(0.35), Da(0.35), n,„h\l3, nbh\ns,a^, k,}. 



Only 



are well constrained us- 
ing SDSS LRGs alone. We marginalize over the other 
parameters, {/3, fJ^/i^, Ws, (t„, fc*}, with the flat priors, 
{(0.1, 0.6), (0.01859, 0.02657), (0.865, 1.059), (0, 500)km 
s"\ (0.09,0. 13)ftMpc"^}, where the flat priors of n^/i^ and 
Us are centered on the measurements f rom WMAP7 and h as 
width of ±7a-wMAP (with awMAP from lKomatsu et al] ( l2010h ). 
These priors are wide enough to ensure that CMB constraints are 
not d ouble counted when our results ar e combined with CMB 
data dChuang. Wang. & Hemanth3 l2012h . We also marginalize 





mean 


a 


lower 


upper 


H{0.35) 


82.1 


5.0 


77.2 


86.9 


Da(0.35) 


1048 


58 


990 


1107 




0.118 


0.017 


0.101 


0.133 




0.0434 


0.0018 


0.0417 


0.0451 


DA{0.35)/rs(zd) 


6.60 


0.26 


6.34 


6.85 


Dvi0.35)/rs(za) 


8.62 


0.25 


8.38 


8.86 


A(0.35) 


0.445 


0.021 


0.425 


0.465 



Table 1. The mean, standard deviation, and the 68% C.L. bounds of 
{H(0.35), 75^(0.35), H^/i^ 7/(0.35) r^^d)^, DA{0.35)/rs(za), 
Dv(0.35)/rsiza), A(0.35)} from SDSS DR7 LRGs. We recommend 
using 7/(0.35) rs{zj^)/c and 7?A(0.35)/rs (z^) for further analysis. The 
unit of 77 is km Mpc"'^. The unit of Da, Dy, and Vg (z^) is Mpc. 



over the amplitude of the galaxy correlation function, effectively 
marginalizing over a linear galaxy bias. 



4 RESULTS 

We now present the model independent measurements of the pa- 
rameters {77(0.35), Da(0.35), f^m^^}, obtained by using the 
method described in previous sections. We also present the de- 
rived parameters including 77(0.35) rs{zd)/c, DA{0.35)/rs{zd), 
Dvi0.35)/r,{za), and yl(0.35), where 



Dv{z) = 



and 



{l + zfDl 



H{z) 



A{z) = Dv{z) 



772 
^0 



(20) 



(21) 



We recommend using {77(0.35) rs{zd)/c, 7)A(0.35)/rs(zd)} in- 
stead of {77(0.35), 7)a(0.35), 0,^11^} because they are more ro- 
bust measurements from this study (see Sec.|5]for more detail). We 
apply our method to the 2D 2PCF of the LasDamas mock cata- 
logs and find that our measurements are consistent with the input 
parameters of the simulations. 



4.1 Constraints on 77(0.35) and Da (0.35) Independent of a 
Dark Energy Model 

Fig. [3] shows one and two-dimensional marginalized contours of 
the parameters, {77(0.35), 7)a(0.35), Qmh^ , 77(0.35) rs(zd)/c, 
7)A(0.35)/r,(zd), Dv{Q.^^)/rs{zd), A(0.35)}, derived in an 
MCMC likelihood analysis from the measured 2D 2PCF of the 
SDSS LRG sample. Table|T]lists the mean, rms variance, and 68% 
confidence level limits of these parameters. Table |2]gives the nor- 
malized covariance matrix for this parameter set. These are inde- 
pendent of a dark energy model, and obtained without assuming a 
flat Universe. 

The constraints on {77(0.35), L>a(0.35), ^^m/l^ 
77(0.35) r4«d)/c, 7)A(0.35)/r,(2d), 7)y(0.35)/r,(zd), 
A(0.35)}, as summarized in Tables [T] and |2l can be used to 
combined with any other cosmological data set to constrain dark 
energy and the cosmological model. We recommend using only 
{77(0.35) r-i,(2;ti)/c, 7)A(0.35)/rs(zd)} since they have tighter 
constraints than {77(0.35), 73a(0. 35)} and are robust in the 
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0.14 
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0.18 
0.16 
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0.12 
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1000 1200 



0.08 0.12 0.16 
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1000 1200 



0.08 0.12 0.16 



0.04 0.045 0.05 





H(0.35) 




1000 1200 

D (0.35) 




0.08 0.12 0.16 




0.04 0.045 0.05 




D (0.35)/r (zj 




D m.35)/r(zj 




Figures. 2D marginalized contours (68% and 95% C.L.) for {/^(O.SS), Da(0.35), Q.mh'^, //(0.35) r4zd)/c, D A{0.?,b) /rs(za), Z)v(0.35)/rs (zd), 
A{0.35)}. The diagonal panels represent the marginalized probabilities. The unit of H is kms~^ Mpc^^. The unit of Dj<^, Dy, and (zj^) is Mpc. 



systematic tests we have carried out (see Sec. |5}- In addition, 
H (0.35) ri,{zd)/c and DA{0.35)/rs{zd) are basically indepen- 
dent to Q,mh? which might not a robust measurement in this study 
(see Sec.[5j. 



10 h ^Mpcx 10 h ^Mpc bin is used if the center of the bin is in 
the scale range of 40/i"^Mpc< s < 120/i^^Mpc. 



The bestfit model from the MCMC likelihood analysis has 
= 112 for 99 bins of data used for a set of 9 parameters 
(including the overall amplitude of the correlation function) and 
the per degree of freedom (x^/d.o.f.) is 1.24. Note that a 



4.2 Validation Using Mock Catalogs 

In order to validate our method, we have applied it to 80 2D 
2PCFs from 80 LasDamas mock catalogs (which are indexed with 
01a-40a and 01b-40b). Again, we apply the flat and wide priors 
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i?(0.35) 


Da(0.35) 




H{0.35)r,{za)/c 


DAiO.35) /r,{za) 


Dv (0.35) /r,{za) 


A(0.35) 


H(0.35) 


1 


-0.4809 


0.7088 


0.7297 


0.0827 


-0.2631 


0.2618 


Da (0.35) 


-0.4809 


1 


-0.6339 


-0.0065 


0.6730 


0.6167 


0.0379 




0.7088 


-0.6339 


1 


0.0867 


0.0888 


0.0427 


0.7042 


//(0.35) rs{zj)/c 


0.7297 


-0.0065 


0.0867 


1 


0.0604 


-0.4104 


-0.1934 


DA(0.35)/r4^d) 


0.0827 


0.6730 


0.0888 


0.0604 


1 


0.8851 


0.6447 


Dv{0.?,b)/rs{zi) 


-0.2631 


0.6167 


0.0427 


-0.4104 


0.8851 


1 


0.6807 


^(0.35) 


0.2618 


0.0379 


0.7042 


-0.1934 


0.6447 


0.6807 


1 



Table 2. Normalized covariance matrix of the measured and derived parameters, {/f(0.35), Z)yi(0.35), f2m/i^, H(0.35) rs{zd)/c, DA(0.35)/rs (z^). 
Dv(0.35)/rs(za),A(0.35)}. 





mean 


a 


input value 


H{0.35) 


81.1 


4.1 


81.79 


Da{0.35) 


1009 


56 


1032.8 




0.121 


0.013 


0.1225 


H(0.35)r,(2d)/c 


0.0434 


0.0020 


0.0434 


DA(0.3b)/rs{za) 


6.26 


0.30 


6.48 


Dv(0.35)/r,(2d) 


8.33 


0.31 


8.51 


A(0.35) 


0.440 


0.019 


0.452 



Table 3. The mean, standard deviation, and the 68% C.L. bounds of the 
distributions of the measured values of {i?(0.35), Da(0.35), Q,Tnh?, 
H(0.35)rsiza)/c, DA(0.35)/r,{2d), (0.35)/r,(2d), A{0.35)} 
from the 2D 2PCF of each of 80 LasDamas mock catalogs (which are in- 
dexed with 01a-40a and 01b-40b). Our measurements are consistent with 
the input values within la, where each a is computed from the 80 means 
measured from the 80 mock catalogs. The unit of H is kms~^ Mpc~^. 
The unit of Da_, Dy, and ^^(z^) is Mpc. 



(±7awMAP7) on Qth^ and Us, centered on the input values of the 
simulation (Qth^ = 0.0196 and — 1). 

Table [5] shows the means and standard deviations of the 
distributions of our measurements of {//(0.35), Da (0.35), 
H{0.35) rs{zd)/c, DA(0.35)/r-,(zd), Dv(0.35)/r-,(zd), 
A(0.35)} from each of the LasDamas mock catalogs (80 total) of 
the SDSS LRG sample. These are consistent with the input param- 
eters, establishing the validity of our method. We also show the 
measurements of H{0.35) rs{zd) /c and DA{0.35)/rs{zd) of each 
mock catalog in Fig.|4]and Fig.|5] One can see the measurement of 
H{0.35) rs{zd)/c and Da{0-S5) /rs{zd) are consistent with the 
input parameters of the simulations. 



4.3 Cross-check with Measurements from Multipoles of the 
Correlation Function 

As a cross-check of our results, we now present the measurements 
from the monopole-quadrupole of the correlation function for com- 
parason. The detail of the method is described in Appendix |B] Ta- 
ble|4]lists the mean, rms variance, and 68% confidence level limits 
of these parameters. The measurements are consistent with those 
from our main method (see Table [T). However, the constraints are 
much weaker (> 8%, which is twice as large as our main results). 
This is most likely due to the fact that the information used in 
the monopole-quadrupole method is much less than what we use 
in our main method (as presented in this paper). It is possible to 
obtain better measurements using multipole method by including 



30a 40a 10b 

index of mock catalogs 



20b 30b 40b 



Figure 4. Measurements of the means of H{0.35) rs(zd)/c from 80 in- 
dividual mock catalogs (indexed as Ola to 40a and 01b to 40b). The black 
solid line shows the mean of these 80 measurements and the blue dashed 
Unes show the range of icr. The red dotted line shows the theoretical value 
computed with the input parameters of the simulations. 




Figure 5. Measurements of the means of DA{0.35)/rs{zd) from 80 in- 
dividual mock catalogs (indexed as Ola to 40a and 01b to 40b). The black 
solid line shows the mean of these 80 measurements and the blue dashed 
lines show the range of ±cr. The red dotted line shows the theoretical value 
computed with the input parameters of the simulations. 
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mean 


(J 


lower 


upper 


H(0.35) 


79.6 


8.8 


70.9 


87.8 




1060 


92 


970 


1150 




0.103 


0.015 


0.088 


0.118 


//(0.35) r,(zd)/c 
r'A(0.35)/r4^d) 


0.0435 
6.44 


0.0045 
0.51 


0.0391 
5.99 


0.0477 
6.90 



Table 4. The mean, .standard deviation, and the 68% C.L. bounds of 
{//(0.35), Da (0.35). C^/i^ H{Q.'ib)rs(za)/c, DA(0.'i5)/rs(zi)} 
from SDSS DR7 LRGs using monopole-quadrapole of the coiTelation func- 
tion. The bin size is 5 h^^Mpc and the scale range is 40 < s < 120 
h-^Mpc. x^/d.o.f is 1.23. The unit of H is kms~i Mpc"^. The unit of 
Da, Dv, and rs(2^) is Mpc. 

hi gher order multipoles of the correlation function; this is explored 
in lchuang & Wand l l2012h . 



5 SYSTEMATIC TESTS 

Table. |5] shows the systematic tests that we have carried out vary- 
ing key assumptions made in our analysis. These include the range 
of scales used to calculate the coiTelation function, the nonlinear 
damping factor, the bin size, and an overall shift in the measured 
correlation function due to a systematic error. 

First, we fix the nonlinear damping factor, fc* — 0.11/iMpc^^, 
and find the results are basically the same. To speed up the compu- 
tation, we fix fc* for the rest of the tests. 

In this study, we marginalize over /3 with a wide flat prior 
(0.1 to 0.6) since our method is not sensitive to /3. We test fix- 
ing the value of /3 to 0.35, which is close to the measurement 
from previous work with s imilar data but using different method 
jCabre & Gaztanagal l2008h . and find that our measurements of 
//(0.35) rs{zd)/c and DAiO-S5) /rs{zd) change by less than 1% 
compared to that of marginalizing over /?. 

We vary the range of the scale and find that (0.35) rs{zd)/c 
and Da (0.35) /r^ (zd) are insensitive to it. However, Qmh^ is sen- 
sitive to the minimum scale chosen which could imply that the 
scale dependent bias or redshift distortion is distorting larger scale 
than we have expected. Therefore, we do not recommend to use 
Qmh^ from this study. In the case of s = 40 — 130/i~^Mpc, 
Da{0-35) /rs{zd) is different from the fiducial result with about 
2(7, which is likely due to systematic errors responsible for the 
anomalously high tail in the spherically-averaged correlation func- 
tion (see, e.g.. lChuang. Wang. & Hemanthj ( l2012h ) on large scales. 

We vary the bin size to 8/i~^Mpcx 8/i~^Mpc and find 
X^/d.o.f.= 1.72, which can be explained by the increase in the 
noise level with the increased number of bins. The number of the 
mock catalogs used to construct the covariance matrix is 160 and 
the number of bins used with bin size = 8/i^^Mpcx8/i~^Mpc is 
159. One can expect the covariance matrix would be too noisy to 
give reasonalbe results. 

We also show that the results are insensitive to the constant 
shift by lowering down the data points of the observed correlation 
function by 0.001 and 0.002. 



6 CONCLUSION AND DISCUSSION 

We have obtained the first measurements of H{z) and 
Da{z) from galaxy clustering data in an MCMC likeli- 



hood analysis. Our constraints for the measured and derived 
parameters, {H{0.35), Da{0.S5), 1^(0.35) r,(zd)/c, 

I>A(0.35)/r,(2d), Dv{0.35)/rs{zd), ^(0.35)}, from the 2D 
2PCF of the sample of SDSS DR7 LRGs are summarized by Ta- 
bles [T] and |2] Our results are robust and independent of a dark en- 
ergy model, and obtained without assuming a flat Universe, and 
represent the first measurements of H{z) and Da{z) from galaxy 
clustering data. 

Our galaxy clustering measurements of H{0.35) rs{zd) /c and 
Da (0.35) /vs (zd) (see Tables[T]and|2j can be used to combine with 
CMB and other c osmological data sets to probe dark energy. In a 
companion paper ( IWang. Chuang. & Mukherieell2012l) . we explore 
the implications of our results for dark energy constraints. 

We recommend using H (0.35) rs{zd)/c and 
DA{0.35)/rs{zd) measured from the SDSS LRGs for com- 
bination with other data sets, since they are tight constraints (4%) 
that are nearly uncorrelated, and robust with respect to tests of sys- 
tematic uncertainties. This is as expected, since (0.35) rs{zd)/c 
and DA{Q-'i^)/rs{zd) correspond to the preferential redshift 
separation along the line of sight, and the preferential angular 
separation in the transverse direction respectively; these in turn 
arise from the BAO in the radial and transverse directions. The 
measurable preferential redshift and angular separations should be 
uncorrelated since they are independent degrees of freedom. On 
the other hand, the measurements of //(0.35) and Da(0.35) are 
mainly determined by the geometrical distortion (i.e. the 2D corre- 
lation function is supposed to be isotropic without considering the 
redshift distortion), so that they are highly correlated (correlation 
coefficient r ~ —0.5). The presence of the BAO (although only 
marginally visible in Fig.l) leads to tight and robust constraints 
on H{0.?,Zi) rs{zd) / c and DA{Q-'i^)/ra{zd)- Since most of the 
constraining power in our analysis comes from fitting the overall 
shape of the galaxy correlation function on quasi-linear scales, and 
not from fitting the BAO peaks, we refer to our measurements as 
galaxy clustering measurements. 

We have validated our method by applying it to the 2D 2PCF 
of the mock catalogs from LasDamas, and finding consistency be- 
tween our measurements and the input parameters of the LasDamas 
simulations for samples (see Table[3}. 

As a cross-check of our results, we have measured H{z) and 
Da{z) using monopole-quadrupole method and find that the re- 
sults are consisent with our main method in this study but the con- 
straints are much weaker. The reason is most likely that the infor- 
mation used by the monopole-quadrupole method is much less than 
what we use in our main method. However, it is still possible to 
improve the constraints by including higher order multipoles. We 
explore this issue in Chuang & Wang (2012). 

Our work has significant implications for future surveys in 
establishing the feasibility of measuring both H{z) and Da{z) 
from galaxy clustering data. In future work, we will optimize our 
method, and apply it to new observational data as they become 
available, and to simulated data of planned surveys to derive robust 
forecasts for dark energy constraints. 
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H{0.35) 



Da(0.35) 



ff(0.35) rs(i .) 



-Pv(0.35) 



A(0.35) 



X /d.o.f. 



fiducial result 
kt = 0.11 
/3 = 0.35 
30 < 3 < 120, fc* = 0.11 
50 < s < 120, fc» = 0.11 
40 < s < 110, fc» = 0.11 
40 < a < 130, fc* = 0.11 
bin size = Sh^^MpcxHh^^Mpc 
shift= 0.001, ft* = 0.11 
shift= 0.002, kj, = 0.11 



,+4.S 
-4.9 
7+4.9 
-5.0 
83.9 ± 5.6 

-4.7 
0+5.5 
-5.4 
80.6 ± 5.1 
a+6.4 

-6.3 
87.9^ 



82.12 
81.7+ 



83. 4^ 
83.9+ 



84.8^ 



83.0 



+ 5.6 
6.0 
+ 5.4 
5.3 
, + 5.4 



1048+J^ 

1051 ± 59 

1008 ± 53 

1038+^? 
-51 

1012+ " 
1087 



9 + 63 
-64 
+ 59 
-60 
g„7+61 
'^°'-60 
1037 ± 60 

1024l:«| 



0.118 ± 0.016 

" "-0.017 
0.132 ± 0.020 

120+°"" 

134+°''2* 
-0.023 

0.115 ± 0.016 
0.115 ± 0.016 

"^^"-0.018 
124+°°1° 

las+o-"^! 



0.0434 ± 0.0017 

0.0434 ± 0.0017 

0430+° ™^'^ 

0437+'' '"'l'' 
"■'J*'^' _0.0018 

0.0428 ± 0.0019 

0.0429 ± 0.0019 

0.0451 ± 0.0026 

0447+"-'"'^^ 
-0.0024 

0433+° °°" 
0-"43J_p QQjg 

0.0435 ± 0.0019 



6.60 ± 0.25 

6.59 ± 0.25 

6.54 ± 0.24 

6 59+°-^^ 
°'°^-0.22 

6.59 ± 0.26 

f. _„ + 0.27 

6 17+°'" 
-0.26 

6.78 ± 0.27 
6.63 ± 0.26 
6.67 ± 0.28 



8.62 ± 0.24 

8.62 ± 0.24 

8.60 ± 0.24 

8 59+°-2'' 
°-^"-0.23 

8 65+°-^'' 
''■'"'-0.23 

8.81 ± 0.27 

8.14 ± 0.28 

„ 7„+0.27 

8-™~0.29 

8 65+°-2= 
"■""-0.26 

8.67 ± 0.28 



0.445 ± 0.020 

0.443 ± 0.020 

0.456 ± 0.021 

0.446 ± 0.016 

0.460 ± 0.024 

454+° °22 
" "^ -0.021 

0.418 ± 0.019 

470+° °" 
" "-0.020 

0.453 ± 0.021 
0.463 ± 0.022 



1.24 
1.24 
1.25 
1.24 
1.06 
1.09 
1.31 
1.72 
1.25 
1.24 



Table 5. This table shows the systematic tests with the damping factor, the scale range, the bin size, and the assumed constant shift from a system- 
atic error (^obsC*) = $true(s)+shift). The fiducial results are obtained by considering the scale range (40 < s < 120 h^^Mpc), the bin size = 
10/i~^Mpcx 10/i~^Mpc, and the damping factor.fc*, marginalized over with the a flat prior (0.09 < fc* < 0.13 hMpc~^). The other results are calcu- 
lated with only specified quantities different from the fiducial one. The unit of H is kms~^ Mpc^^. The unit of Da, Dy, and rs(z^) is Mpc. The unit of 
fc* is hMpc~^. 
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APPENDIX A: ALGORITHM OF SMOOTHING THE 
COVARIANCE MATRIX 

The original covariance matrix, Ci^ , is computed by Eq. l |13t . Since 
the correlation function is in two dimension, we could re-label 
the covariance matrix as a function of the indexes of two bins, 
C{(Ji,Tii,(jj,-nj). The covariance matrix would be noisy since it 
is constructed with a small number of the mock catalogs compar- 
ing to the number of bins used. One might not obtain converging 
results while applying MCMC analysis on the observed data or an 
individual mock catalog, especially considering larger scale range 
or larger number of bins. We find that smoothing the covariance 
matrix could solve the problem. We have also checked that the 
smoothing procedure would not introduce bias by comparing the 
results from applying the smoothed and the original covariance ma- 
trix on the averaged correlation function from 160 mock catalogs. 
The concept of our smoothing procedure is that the new value of 
an element of the covariance matrix would be determined by its 
original value and the values of its neighbor elements. For exam- 
ple, to smooth an array, f[ri\, we assign the new value at the index, 
n', by /[n'] = (1 - p) ■ f[n'] + p ■ (/[n' - 1] + f[n' + l])/2 
where ^ p ^ 1 . The goal is to make the value be closer to the 
mean of the neighbors. Notice that while one of the neighbors is 
not available (i.e. f[n'] is the first or last element of the array), we 
let f[n'] be fixed since the mean of the neighbors is not available. 
The algorithm can be expressed by eq. JAlt 

While eq. JAlt can be applied on most of the elements of the 
covariance matrix, there are some special cases as decrible below. 

(1) Diagonal elements: these elements are only determined by the 
nearby diagonal elements by eq. ( IA2t . since the diagonal elements 
are supposed to be larger than their off-diagonal neighbors and 
should not be smoothed using the latter. 

(2) First off-diagonal elements: these elements are only determined 
by the first off-diagonal elements nearby, i.e. C((Ji+As, tt;, tJi, tt^) 
would be assigned a new value by eq. iA3t , since a first off- 
diagonal element could be very different from its neighboring di- 
agonal elements and should not be smoothed using the latter. 



APPENDIX B: MEASURING H AND Da WITH 
MULTIPOLES OF THE CORRELATION FUNCTION 

Using monopole-quadrupole of power spectrum to break 
the degeneracy of H (z ) an d Da(z) is i ntrodu ced by 
IPadmanabhan a nd White! ( 120081) . iKazin et alj j201ll) tested 
the method of monopole-quadrupol e of correlation fun ction with 
mock catalogs. The method used in iKazin et alj i201lh is similar 
but not exactly the same as our method. We describe our method 
below. 

First, we compute the 2D correlation function with bin size, 
1 fe~^Mpcx 1 /i^^Mpc. Then, we compute the monopole and 
quadrupole by 

?o(*) = -r; — ; tt- r- — ; (Bl) 

Number oi bms used m the numerator 

and 

Number of bins used in the numerator ' 

where 
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, TTj) = (1 — p) ■ C(cri,1Ti,Crj,TTj) + 



(C{cri + As, TTi, Oj , TTj) + C((Ti — As, TTi , (Tj , vTj)), if both elements left are available 
{C'((Ti,7Ti + As, Oj , TTj) + C{cri,iTi — As, (7j,TTj)), if both elements left are available 
{C'{cri , TTi, Oj + As, ■Kj) + C((Ti , TTi, Oj — As, TTj)), if both elements left are available 
(C(o"i , TTi, (Tj , TTj + As) + C{cTi , TTi, cTj , TTj — As)), if both elements left are available 



where m is the number of the neighbor elements used which should be 0, 2, 4, 6, or 8, and As is the size of the bin in one direction. We use p : 
As = 10 Mpc/h in this study, but p = if m = 0. And then, we iterate this procedure for 10 times. 



: 0.01 and 



C(CTi, TTi, (Ti, TTi) 



(1 — p) ■ C{c!i,-Ki,ai,-Ki) 



(A2) 



P_ 

m 



(C(<7„ 



■ As, TTi, 

H + As, 



■ As,7rj) + C{tJi 

li + As) + C{CJ^, 



As, TTi , (Ti — As, TTi)), if both elements left are available 
'i — As, cFi,iTi — As)), if both elements left are available 




Figure Bl. Measurement of monopole of the correlation function from 
SDSS DR7 LRGs in a redshift range 0.16 < 2 < 0.44. The enor bars 
are the square roots of the diagonal elements of the covariance matrix we 
have derived from mock catalogs. 



Figure B2. Measurement of quadrupole of the correlation function from 
SDSS DR7 LRGs in a redshift range 0.16 < ^ < 0.44. The eiTor bars are 
the square roots of the diagonal elements of the covariance matrix we have 
derived from mock catalogs. 



and As = 5 h^^Mpc and the scale range used is 40 < s < 120 
h~^Mpc in this study. Therefore, there are 32 data points for 
monopole-quadrupole method here, with 16 measurements each for 
monopole and quadrupole respectively. Fig. IB II and IB2l show the 
measurements of the monopole and quadrupole of the correlation 
function from the observed galaxy sample. 

Just like our main method, the covariance matrix (for the 32 
monopole and quadrupole measurements) is constructed from the 
mock catalogs. The theoretical multipoles are computed using Eqs. 
dBTT l and ( |B2] |. 

Now, by exploring the same parameters and ranges using 
MCMC (with given by Eq.fTel). one could measure H[z) and 
Da{z) following the same steps as our main method in this study. 
The results are shown in sec. 14.31 (see Table 4). Like our main 
method, our monopole-quadrupole method has only one approx- 
imation, that the observed 2D correlation function using differ- 
ent fiducial models can be converted from one to another with 
two stretching fac tors. However, the multipole method tested in 
lKazinetalJj201ll) neglected some additional terms while measur- 
ing the stretching factors. Although the effect of the neglected terms 
might be small, it could be completely avoided by using the method 
described here. 
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(C(<Ti + 2As,7ri,ffi 
- AsjTTi + As,ai,iri ■ 



(A3) 



- As, TTj) + C((Ti , TT j , (Tj — As, TTj)), if both elements left are available 

■ As) + C((Ti +As,7ri — As,(Ti,7ri — As)), if both elements left are available 
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